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METHOD AND APPARATUS FOR 
EVALUATING STORED CHARGE IN AN 
ELECTROCHEMICAL CELL OR BATTERY 

BACKGROUND OF THE INVENTION 
5 The present invention relates to testing of 

storage batteries. More specifically, the invention 
relates to evaluating stored charge in an 
electrochemical cell or battery. 

Stored charge is an important parameter in 

10 many applications of electrochemical cells and 
batteries. With traction batteries, stored charge 
represents an electric vehicle's fuel supply and thus 
determines how far the vehicle can travel before 
recharging. With stationary standby batteries, the 

15 level of stored charge determines how long a critical 
load can continue to function in the event of a power 
failure or disconnection from the ac mains. In 
automotive applications, stored charge determines the 
length of time that the lights and accessories can be 

20 operated when the - engine is off, or when the charging 
system has malfunctioned. 

With lead-acid batteries, relative stored 
charge, or state-of -charge (SOC) , has been 
traditionally evaluated by observing either the 

25 . battery's open-circuit voltage, or the specific gravity 
of the battery's electrolyte. However, neither of 
these measurements yields an absolute determination of 
the amount of stored charge. Furthermore, specific 
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gravity measurements are messy and altogether 
impossible to perform on sealed lead-acid cells; and 
open-circuit voltage is difficult to determine under 
load conditions, and is imprecisely related to SOC 
5 since it is greatly affected by both "surface charge" 
and temperature. 

Because of these problems, several 
techniques for correcting the voltage of lead-acid 
batteries to obtain SOC have been proposed. These 

10 include the techniques described by Christianson et al . 
in U.S. Patent 3,946,299, by Reni et al . in U.S. Patent 
5,352,968, and by Hirzel in U.S. Patent 5,381,096. 
However, such voltage correction methods are not very 
accurate. Furthermore, with electrochemical systems 

15 other than lead-acid, they may be of little help since 
battery voltage often bears very little relationship to 
stored charge. 

Because of the problems with traditional 
methods for determining relative charge (SOC) , many 

20 techniques based upon measuring ac impedance have been 
suggested. For example, U.S. Patent 3,984,762 to 
Dowgiallo purports to determine SOC directly from the 
phase angle of the complex impedance at a single 
frequency. In U.S. Patent 4,743,855, Randin et al . 

25 assert that SOC can be determined "from the argument . 
(i.e., phase angle) of ' the difference between complex 
impedances measured at two different frequencies. 
Bounaga, in U.S. Patent 5,650,93 7, reportedly 
determines SOC from measurements of the imaginary part 
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of the complex impedance at . a single frequency. 
Finally, " Basell . et al. f in U.S.- Patent 5,717,33-6 
purport to determine SOC from the rate of change of 
impedance with frequency at low frequency. However, the 
fact that none of these ac impedance methods has gained 
wide acceptance suggests that they may not be 
altogether satisfactory methods for determining SOC. 

The 'absolute stored charge or amp -hour 
capacity of batteries has been traditionally measured 
by timed-discharge tests. However, because of the 
expense and the time involved in ^performing such tests, 
ac techniques for determining amp-hour capacity have 
been proposed. Sharaf, in U. S. Patent 3,808,522, 
teaches a method for determining the capacity of a 
lead-acid battery from measurements of its ac internal 
resistance. Yang, in U.S. Patent 5,126,675, also uses 
measurements of internal resistance to determine 
battery capacity. Muramatsu reports, in U.S. Patent 
4,678,998, that he can determine both the remaining 
amp-hour capacity and the remaining service life of a 
battery from measurements of the ac impedance magnitude 
at two different frequencies. Fang, in U.S. Patent 
5,241,275, teaches a method for ' determining remaining 
capacity from complex impedance measured at two or 
three frequencies in the range from 0.001 to 1.0 Hz. 
Finally, Champlin, in U.S. Patent 5,140,269, has shown 
that percent capacity can be determined from the 
measured dynamic conductance at a single frequency if 
the dynamic conductance of a reference, fully- charged, 
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identically constructed, new battery is known. This 
method, although . quite accurate, requires that aprioi 
data be available. 

SUMMARY OF THE INVENTION 
A testing device applies time -varying 
electrical excitation to a cell or battery and senses 
the resulting time-varying electrical response. 
- Computation circuitry within the device uses voltage 
and current signals derived from the excitation and 
response signals as inputs and computes values of 
elements of an equivalent circuit representation of the 
cell or battery. In one aspect, the relative charge 
(SOC) of the cell or battery is calculated from the 
value of the conductance component G of a particular 
parallel G-C subcircuit of the equivalent ci.rcuit . In 
another, the absolute charge (Ah) contained in the cell 
or. battery is calculated from the value of the 
capacitance component C of a different parallel G-C 
subcircuit. In other aspects, relative or absolute 
charge values are then either displayed to the user or 
are used to control an external process such as 
charging of the battery. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of a device for 
evaluating stored charge in an electrochemical cell or 
battery according to the present invention. 

FIG. 2 depicts a six-element small signal 
equivalent circuit representation ' of a particular 
automotive storage battery. 
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FIG. 3 is a plot of the variation of the 
three subcircuit time-constants defined in FIG. 2 as 
functions of the charge removed from the battery. 

FIG. 4 is a plot of measured and theoretical 
values of conductance G3 defined in FIG. 2 as functions 
of the charge removed from the battery. 

FIG. 5 is a plot of the conductance, or real 
part of the admittance, of an automotive battery 
measured at both 22 Hz and 8 Hz as functions of the 
mean specific gravity of the battery. 

FIG. G is plot of measured values of 
capacitance C2 defined in FIG. 2 as a function of the 
charge removed from the battery. 

FIG. 7 is a plot showing experimental and 
theoretical values of the amp-hour capacity of six 
fully-charged batteries as functions of the value of 
capacitance C2 defined in FIG. 2. 

FIG. 8 is a circuit representation of the 
parallel G3-C3 subcircuit showing its admittance Y3 . 

FIG. 9 is a circuit representation of the 
parallel G2-C2 subcircuit showing its admittance Y2 . 
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A method and apparatus for quickly and 
accurately determining both relative charge (SOC) and 
absolute charge (Ah) that is unaffected by "surface 
charge" and does not require apriori data would be of 
great value. The present invention addresses this 
need. It is based upon teachings disclosed in pending 
U.S. Patent Application Ser. No. 09/152,219, filed 
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September 11, 1998 and entitled "METHOD AND APPARATUS 
FOR MEASURING COMPLEX IMPEDANCE OF CELLS AND 
BATTERIES" and pending U.S. Patent Application 
09/151,324, filed September 11, 1998, entitled 
5 "METHOD AND APPARATUS ' FOR DETERMINING BATTERY 
PROPERTIES FROM COMPLEX IMPEDANCE ADMITTANCE" which 
are incorporated herein by referenced 

Figure 1 discloses a block diagram of 
apparatus for evaluating stored charge according to 

10 the present invention. Apparatus of this type is 
fully disclosed in pending U.S. Patent Application 
Ser. No. 09/152,219 and pending U.S. Patent 
Application 09/151,324. Measuring circuitry 10 

electrically couples to cell/battery 20 by means of 

15 current -carrying contacts A and B and voltage-sensing 
contacts C and D. Measuring circuitry 10 passes a 
periodic time-varying current i(t) through contacts A 
and B and senses a periodic time-varying voltage v(t) 
across contacts C and D. By appropriately processing 

20 and combining i(t) and v.(t) , measuring circuitry 10 
determines real and imaginary parts of a complex 
' parameter, either impedance Z or admittance Y, at a 
measuring frequency f k ; where f k is a discrete 
frequency contained in the periodic waveforms of both 

25 i (t) and v(t) . 

Control circuitry 3 0 couples to measuring 
circuitry 10 via command path 40 and commands 
measuring circuitry 10 to determine the complex 
parameter of cell/battery 20 at each one of n 
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discrete measuring frequencies, where n is an 'integer 
number. This action defines 3n experimental 

quantities: the values of the n measuring frequencies 
and the values of the n imaginary parts and n real 
parts of the complex parameter at the n measuring 
frequencies . 

Computation circuitry 50 couples to 
measuring circuitry 10 and to control circuitry 30 
via data paths SO and 70, respectively, and accepts 
the 2n experimental values from measuring circuitry 
10 and the values of the n measuring frequencies from 
control circuitry 30. Upon a "Begin Computation" 
command from control circuitry 3 0 via command path 
80, computation circuitry 50 uses algorithms 
disclosed in D. S. Patent Application 09/151,324 to 
combine these 3n quantities numerically to evaluate 
2n elements of an equivalent circuit representation 
of the cell/ battery. Computation circuitry 50 then 
calculates the relative charge (SOC) and/or the 
absolute charge (Ah) of the cell/battery from values 
of .particular elements of this circuit 
representation. Finally, computation circuitry 50 
outputs the computed result to the user on display 90 
and/or uses the result to control a process 100 such „ 
as a battery charger. 

In practice, a microprocessor or 
microcontroller running an appropriate software 
program can perform the functions of both control 
circuitry 30 and computation circuitry 50. 
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Figure 2 discloses a six-element equivalent 
circuit representation of a typical automotive 
storage battery. This circuit representation was 
evaluated using apparatus of the type disclosed in 
5 FIG 1 with n=3 by employing algorithms disclosed in 
U.S. Patent Application 09/151,324. The three 

measurement frequencies were 5 Hz, 7 0 "Hz, and 10 0 0 
Hz. One notes that the' n=3 equivalent circuit 
comprises three subcircuits: 
10 • A series Gl-Ll subcircuit. 

• A parallel G2-C2 subcircuit. 

• A parallel G3-C3 subcircuit. 

One notes further that the three subcircuits 
are characterized by having very different time 
15 constants. The shortest time constant, 

tj =L1 Gl = 93.5jaS , belongs to the series Gl-Ll 
subcircuit. The next longest time constant, 

i 2 = C2/G2 = 2.22raS , belongs to the parallel G2-C2 
subcircuit; and the longest time-constant, 

20 x 3 =C3/G3 = 41.6mS , belongs to the parallel G3-C3 
subcircuit. Accordingly, the three subcircuits 

represent quite different physical processes and can be 
differentiated from one another by their time 
constants . 

25 Figure 3 is a logarithmic plot of the three 

time constants defined above as functions- of charge 
(ampere-hours) removed from the battery. One notes 
that the three time constants remain widely separated 



WO 01/16615 



PCTVUS00/40816 



as charge is removed, and that the longest of the 
three, t 3/ is nearly independent of state -of -charge . 

Figure 4 discloses the variation of 
conductance G3 with charge (amp-hours) removed from the 
5 battery. One sees that G3 approaches minimum near full 
charge and again near full discharge, while reaching 
maximum at about 50% state-of -charge . This variation 
is consistent with a theoretical model that associates 
the G3-C3 subcircuit with a linearized, small-signal, 
10 representation of the nonlinear electrochemical 
reaction occurring at the negative plates . For such a 
model, conductance G3 would be proportional to both the 
number of reaction sites available for the charge 
reaction (i.e., PBSQ . sites) and also the number of' 
15 sites available for the discharge reaction (i.e., Pb 
sites) . Accordingly, G3 would be proportional to the 

product (SOC)-(l-SOC) . Figure 4 shows a theoretical 
curve based upon this assumption. The agreement 
between the theoretical curve and the experimental 

20 points is seen to be excellent. 

The observed variation of G3 with SOC can be 
exploited to determine SOC from measurements of G3 . 
Inverting the . theoretical G3 (SOC) curve leads to a 
quadratic equation for SOC(G3). This • quadratic 

25 equation has two roots, identified as SOC + and SOC" for 
G3=400S in FIG. 4. The SOC* root in FIG 4 corresponds 
to approximately 85% SOC and the SOC" root corresponds 
to about 15% SOC. Thus, if one knows which root is the 
correct root, one can determine SOC from measurements 
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of G3 . Introducing a single piece of auxiliary 
information such as battery voltage, single- frequency 
conductance, or the value of one of the other elements 
of the equivalent circuit can readily ■ identify the 
5 correct root . 

Figure 5 discloses plots of measurements of 
conductance, i.e. the real part of the complex 
admittance, obtained at two different frequencies 8 
Hz and 22 Hz. These results were obtained by 

10 sequentially discharging a battery while measuring the 
conductance and the specific gravity after each 
discharge. The data are plotted as functions of the 
battery electrolyte's mean specific gravity. Although 
the effect is not as pronounced, one notes the same 

15 type of behavior at both frequencies as was noted above 
with reference to FIG. 4. That is, the measured 
conductance is smallest at both full charge and at full 
discharge, and reaches maximum' near 50% SOC. This 
suggests that if accuracy is not too important, one may 

20 be able to determine an approximation to SOC from a 
simple analysis of measurements of conductance alone 
obtained at one or more appropriately chosen 
frequencies. Such an analysis would be much simpler 
than the rigorous mathematical analysis disclosed in 

25 pending U.S. Patent Application 09/151,324 and would 
take advantage of the fact, disclosed above in FIG. 3, 
that the time constants of the subcircuits are widely 
separated so that the subcircuits have little 
interaction. 



WO 01/16615 



PCT/USOO/40816 



-11- 

Figure 6 discloses a plot of measured values 
of capacitance C2 as a function of charge (amp-hours) 
removed from the battery. Note that this capacitance 
is observed to decrease monotonically with charge 
5 removed. I have found the value of capacitance C2 to 
be an excellent indicator of the absolute charge (Ah) 
contained in the battery. As such, its measurement 
could be used to implement a battery "fuel gauge' 7 . To 
illustrate this possibility, "EMPTY and "FULL" 

10 indications have been added to the plot of FIG. 6. 
Although the relationship between C2 and remaining 
charge is not exactly a linear relationship, it appears 
to be fairly close to linear. One notes that, after 
all, many automobile fuel ' gauges are not exactly linear 

15 either. 

Figure 7 discloses the total absolute charge 
in amp-hours as a function of measured capacitance C2 
for six fully-charged batteries of varying sizes. The 
experimental amp-hour points plotted for each battery 

20 were obtained from actual timed discharge tests 
obtained with a constant discharge current of 25 
amperes, to a final terminal voltage of 10.5 volts. 
Also plotted in FIG. 7 is an experimental linear 
equation of the form Ah = 22 + 8.2- C2 . One notes chat the 

25 experimental points agree with the linear equation 
within approximately ±3 ampere-hours. This 
interesting result implies that the absolute charge 
capacity in ampere-hours can be quite accurately 
determined . from measurements of C2 without the 
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necessity of performing- costly and time-consuming 
timed-discharge -tests. 

Figures 8 and 9 illustrate the aspects of the 
invention. Figure 8 illustrates the G3-C3 subcircuit 
5 and shows that the complex admittance of this parallel 
subcircuit is Y3 = G3-H j©C3 . Thus, my discussion above 
actually discloses a relationship existing between the 
real part of Y3 and the SOC of the battery. Similarly, 
FIG. 9 illustrates the G2-C2 subcircuit and shows that 
10 its complex admittance is Y2 = G2 + jcoC2 . Thus, my 
discussion above actually discloses a relationship 
existing between absolute charge Ah)- and the imaginary 
part of Y2 at a particular frequency co . Although it 
is true that complex Z and complex Y are reciprocals of 
15 one another, the same is not true of their real parts 
(R,G) or their imaginary parts K,B) . Accordingly, the 
results of any ac measurement must be expressed in 
complex admittance form - not complex impedance form - 
in order to observe the important relationships that I 
20 have disclosed herein. 

This fundamental result sets my " work apart 
from the prior-art work of those sited above. Previous 
workers Dowgiallo (in U.S. Patent 3,984,762), Randin et 
al. (in U.S. Patent 4,743,855), Bounaga (in U.S. Patent 
25 5, 650,937), Basell et al . (in U.S. Patent 5,717,336), 
Sharaf (in U.S. Patent 3,808,522), Yang (in U.S. Patent 
5, 126, 675), Muramatsu (in U.S. Patent 4,678,998), and 
Fang (in U.S. Patent 5,241,275) have all attempted to 
derive either SOC or Ah information from the magnitude 
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of the complex impedance, or from components of complex 
impedance measured at one or more frequency. However, 
none of these workers has suggested that the real or 
imaginary parts of a complex admittance were of any 
5 interest at all. 

Although my disclosure has relied upon 
particular apparatus and algorithms previously 
disclosed in U.S. Patent Applications Ser. No. 
09/152,219 and Ser. No. 09/151,324, other methods will 

10 be apparent to one skilled in the arts. For example, 
one can employ bridges or other types of apparatus to 
measure complex admittance or . impedance . Furthermore, 
if accuracy is not a strict requirement, one can take 
advantage of the fact that the various time constants 

15 are widely separated from one another and simply assume 
that the subcircuits are not coupled. Within this 
approximation, C2 and C3 are treated as short circuits 

at frequencies near f 01 =1/2711, , LI and C3 are treated as 

short circuits at frequencies near f o: = 1/2t:t : , and at 

20 frequencies near f 03 = l/2m 3 , LI are treated as a short 
circuit while C2 is treated as an open circuit. Thus, 
with some batteries, ' it is possible to obtain 
satisfactory results with a very simple analysis of 
measurements at two or three frequencies. With certain 

25 batteries, it . is even possible to obtain useful 
approximations to C2 or G3 from measurements of complex 
Y or Z = 1/Y obtained, at a single, appropriately 
chosen, frequency. Workers skilled in the art will 
recognize that these and other variations may be made 
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in form and detail without departing from the true 
spirit and scope of my invention. 
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WHAT IS CLAIMED IS: 
1. Apparatus for determining state-of -charge of an 
electrochemical cell or battery comprising: 

electrical excitation circuitry adapted to 
apply time-varying electrical 

excitation to said cell or battery; 
response sensing circuitry configured to 
sense time-varying electrical response 
of said cell or battery generated as 
a result of said time-varying 
electrical excitation; and ' 
computation circuitry responsive to said 
time-varying excitation and to said 
time -varying response and adapted to 
calculate the value of the real part 
of the complex admittance of a 
particular parallel subcircuit 

comprising part of an electrical 
circuit representation of said cell or 
battery, said real part' characterized 
by varying with state -of -charge and 
reaching maximum value at 

approximately 50% state-of -charge, 
said computation circuitry further 
adapted to evaluate said state-of - 
charge from the computed value of said 
real part of said complex admittance. 
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2 . The apparatus of claim 1 including display 
circuitry and wherein the computation circuitry 
provides the computed value of said state -of - charge 
to the display. 

3 . The apparatus of claim 1 wherein said 
computation circuitry couples to an external process 
device and said external process device is controlled 
by said computation circuitry in accordance with the 
computed value of said state-of -charge . 

4 . The apparatus of claim 1 wherein said 
electrical circuit representation comprises at least 
three subcircuits having differing time constants and 
the time constant of . said particular parallel 
subcircuit is longer than the time constants of at 
least two other subcircuits. 

5. The apparatus of claim 1 wherein said 
electrical circuit representation comprises at least 
two subcircuits having different time constants and 
the time constant of said particular parallel 
subcircuit is longer than the time constant of at 
least one other subcircuit. 

6. The apparatus of. claim 1 wherein said 
computation means employs auxiliary data to 
differentiate between two alternative computed values, 
of said state-of -charge . 

7. -Apparatus for determining the absolute 
charge contained in an electrochemical cell or 
battery comprising: 
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electrical excitation circuitry adapted to 
apply time -varying electrical 

excitation to said cell or battery; 

response sensing circuitry configured to 
sense time-varying electrical response 
of said cell or battery generated as 
a result of said time -varying 
electrical excitation; and 

computation circuitry responsive to said 
time-varying excitation and to said 
time-varying response and adapted to 
compute the value of a capacitive 
element contained in a particular 
parallel subcircuit comprising part of 
an electrical circuit . representation 
of said cell or battery, said 
computation circuitry further adapted 
to evaluate said absolute charge from 
said computed value of said capacitive 
element . 

8. The apparatus of claim 7 including a 
display device coupled to said computation circuitry 
and wherein the computation circuitry provides said 
evaluated absolute charge. 

9. The. apparatus of claim 8 wherein said 
display device comprises a battery fuel gauge. 

10. The apparatus of claim 7 wherein said 
computation circuitry couples to an external process 
device and said external process device is controlled 
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by said computation circuitry in accordance with the 
computed value of said absolute charge. 



electrical circuit representation comprises at least 
three subcircuits having different time constants and 
the time constant of said particular parallel 
subcircuit is neither the longest nor the shortest 
said time constant. 

12 . ■ The apparatus of claim 7 wherein said 
electrical circuit representation comprises at least 
two - subcircuits having different time constants and 
the time constant of said particular parallel 
subcircuit is longer than the time constant of at 
least one other subcircuits. 

13. A method for determining the state-of- 
charge of an electrochemical cell or battery 
comprising : 



11. 



The apparatus of claim 7 wherein said 



applying 



time-varying 



electrical 



excitation 



to 



said 



cell 



or 



battery; 



sensing 



time -varying 



electrical 



response 



to 



said 



electrical 



excitation; 



evaluating the real . part of admittance 
of a particular parallel G-C 
subcircuit of an equivalent 
circuit representation of said 
cell or battery from said time- 
varying electrical excitation and 
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said time -varying electrical 
response ; and, 
determining said state-of -charge from 
said real part of admittance' of 
• said G-C subcircuit. 

14. A method as in claim 13 wherein determining 
the ' state-of -charge from said real part of admittance 
includes using auxiliary data to select one value 
from two alternative values of said state-of -charge 
determined from said real part of admittance. 

15 . A method as in claim 14 including 
displaying said selected one value of said state-of - 
charge . 

16. A method as in claim 14 including 
controlling an external process device in accordance 
with said selected one value of said state-of -charge . 

17. A method for determining the value of the 
absolute charge contained in an electrochemical cell 
or battery comprising: 

applying time-varying electrical 
excitation to said cell or 
battery; 

sensing time-varying electrical 

response to said electrical 
excitation; 

evaluating the capacitance of a 
capacitive element of a 

particular parallel G-C 

subcircuit of an equivalent 
circuit representation of said 
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cell or battery from said time- 
varying electrical excitation and 
. said time-varying electrical 
response; and, 
evaluating said absolute charge from 
the value of said capacitive 
element . 

18. The method of claim .17 including displaying 

■said evaluated absolute charge to a user. 
13- The method of claim 17 including 

controlling an external process device in accordance 
with said evaluated absolute charge. 
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